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ABSTRACT 

We examine the likelihoods of different cosmological models and cluster evolutionary 
histories by comparing semi-analytical predictions of X-ray cluster number counts 
to observational data from the ROSAT satellite. We model cluster abundance as a 
function of mass and redshift using a Press-Schechter distribution, and assume the 
temperature T(M, z) and bolometric luminosity Lx{M, z) scale as power laws in mass 
and epoch, in order to construct expected counts as a function of X-ray flux. The 
L x — M scaling is fixed using the local luminosity function while the degree of evolution 
in the X-ray luminosity with redshift Lx oc(l + z) s is left open, with s an interesting 
free parameter which we investigate. We examine open and flat cosmologies with 
initial, scale-free fluctuation spectra having indices n = 0, — 1 and —2. An independent 
constraint arising from the slope of the luminosity-temperature relation strongly favors 
the n=— 2 spectrum. 

The expected counts demonstrate a strong dependence on Qq and s, with lesser 
dependence on Ao and n. Comparison with the observed counts reveals a "ridge" of 
acceptable models in Qq — s plane, roughly following the relation s ~ 6f2o and spanning 
low-density models with a small degree of evolution to Q = 1 models with strong 
evolution. Models with moderate evolution are revealed to have a strong lower limit 
of f^o f> 0.3, and low-evolution models imply that Qq < 1 at a very high confidence 
level. We suggest observational tests for breaking the degeneracy along this ridge, and 
discuss implications for evolutionary histories of the intracluster medium. 

Key words: cosmology: theory - cosmology: observation - clusters: galaxies: general 
- dark matter 



1 INTRODUCTION 

Rich galaxy clusters are the youngest virialized objects ex- 
tant, and as such, they provide a unique source of informa- 
tion about our universe. Observations of clusters in the X- 
ray band provide useful information about large-scale struc- 
ture and galaxy formation. Detailed X-ray images provide 
information about the structure of individual clusters, and 
surveys such as the ROSAT all-sky survey (RASS) pro- 
vide data on the general cluster population. In particu- 
lar, counts of clusters as a function of their X-ray flux — 
the logN-logS relation — provide an avenue for exploring 
cosmic evolution of clusters. Recently, several groups have 
pushed the logN-logS relation to fluxes nearly an order 



of magnitude fainter than previous determinations (Rosati 
et al.1995; Rosati & Delia Ceca 1997; Jones et al.1997). In 
this paper, we compare these data to predictions of viable 
cosmological and evolutionary models for clusters, with the 
goal of defining the range of combined cosmology and X-ray 
evolution consistent with the current data. Because the data 
underconstrain the theory, a range of degenerate models 
emerge, but these can be distinguished by upcoming ob- 
servational tests of the high redshift cluster population. 

In order to proceed, we assume the distribution of clus- 
ters as a function of mass and redshift is accurately de- 
scribed by a Press-Schecter (1974) abundance function. To 
investigate the behavior of a wide class of structure forma- 
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tion models, we employ scale-free initial Gaussian perturba- 
tions with spectral indices n equal to 0, —1, and —2 in the 
cluster mass regime. (The effect of curvature in the spec- 
trum is discussed below.) We allow fio to vary and inves- 
tigate both open models without a cosmological constant 
(Ao = A/3 H 2 = 0) and flat models where ilo + Xo = 1- We 
assume a value of H = 50 km s _1 Mpc -1 throughout this 
work, in accordance with the traditional treatment of X-ray 
data, but most of our analysis is independent of this choice. 

The Press- Schecter formalism predicts the number den- 
sity of collapsed dark matter halos as a function of their mass 
and redshift. Since what we observe is the band-limited X- 
ray luminosity (Lx(E)), we need to find a way to relate 
these quantities to one another. Rather than adopt a specific 
"microphysical" model to describe the relationship between 
mass and bolometric luminosity, we assume that it can be 
adequately described by a power law in the mass regime of 
interest, and then fit the free parameters of the power law 
using the local X-ray luminosity function (XLF). The emis- 
sivity of the gas is modeled as a thermal Bremsstrahlung 
spectrum, which is integrated to determine the fraction of 
the bolometric luminosity falling in the appropriate energy 
band. We use the ROSAT Brightest Cluster Sample (BCS) 
compiled by Ebeling et al.(f997) to constrain the parame- 
ters of the fit. 

Since X-ray luminosity is proportional to the squared 
density of ions in the intracluster medium (ICM), and since 
the early universe was denser than it is today, it is reasonable 
to expect that distant clusters may be stronger X-ray emit- 
ters despite having lower overall masses. Arguments based 
on self-similarity predict evolution in the bolometric X-ray 
luminosity of the form L x oc M 4/3 (l + z) 7/2 (Kaiser 1986), 
implying strong positive evolution in the luminosity of ob- 
jects at a fixed mass. However, self-similarity is not par- 
ticularly well justified either obervationally or theoretically; 
models which invoke a minimum central entropy in the clus- 
ter gas fare better in many respects (Evrard & Henry 1991; 
Kaiser 1991). 

The issue of evolution in the XLF has been hotly pur- 
sued among observers; whether or not any evolution exists, 
however, still seems an open question. Early work on the Ex- 
tended Medium Sensitivity Survey (EMSS) cluster sample 
by Henry et al. (1992) and Gioia et al. (1990a) found intrigu- 
ing evidence for evolution in the cluster population. These 
results sparked interest but were preliminary analyses that 
suffered from occasional misclassifications and the vagaries 
of small-number statistics. These same clusters were recently 
reanalyzed with new X-ray observations by Nichol et al. 
(1997), who found that the data were consistent with no 
evolution in the cluster population out to redshifts of about 
0.3. The BCS set of 199 clusters mentioned above also shows 
no evidence of evolution in the XLF at low redshifts. More 
distant observations include those of Castander et al. (1995), 
Luppino & Gioia (1995), and Collins et al.(1997). Luppino 
and Gioia have been collecting high-redshift clusters from 
the EMSS sample and finding number densities higher than 
the standard fi = 1 CDM model predicts, which could indi- 
cate either a low-density universe or the existence of a strong 
degree of evolution. Castander et al.have found a dearth of 



high-redshift clusters relative to a simple no-evolution pre- 
diction based on integration of the local XLF. Collins et 
al., however, compiled a sample of high-reshift clusters from 
the SHARC (Serindipitous High-redshift Archival ROSAT 
Cluster) survey, and found significantly higher numbers at 
z > 0.3. 

Given the lack of consensus in this debate, we allow for 
freedom in how luminosity scales with redshift at a fixed 
mass by writing Lx oc (1 + z) s . The cosmological models 
described above are then examined under varying degrees 
of evolution in the X-ray luminosity, yielding four interest- 
ing free parameters for the models considered in this paper: 
fio, Ao, n, and s. Armed with the above assumptions, we 
extrapolate the local cluster abundances to high redshifts 
and predict the surface densities of clusters at low fluxes. 
Comparing predictions to recent observational determina- 
tions of the cluster logN-logS (Rosati & Delia Ceca 1997; 
Jones et al.1997) allows us to make statements about the 
relative probabilities of the models. 

In §2, we describe the mathematical model used and 
the process of fixing a subset of parameters with local ob- 
servations. In §3, we describe our method of predicting the 
logN-logS and the status of current measurements of this 
quantity. We also discuss the likelihoods of individual mod- 
els and the cosmological constraints which can be obtained 
from analysis of the logN-logS alone. In §4, we go on to ex- 
plore further ways of discriminating among cosmologies, by 
making use of the luminosity-temperature relation and red- 
shift distributions of flux limited samples. Finally, in section 
5, we sum up our results and suggest future directions for 
theory and observation in this rich field. 



2 FIXING THE MODEL 
2.1 Theoretical Framework 

The method described here is similar to that put forth in 
Evrard & Henry (1991, hereafter referred to as EH91). The 
first step in predicting the number of observable X-ray clus- 
ters at a given flux limit is to model the distribution of these 
objects. Given that X-ray clusters correspond to virialized 
dark matter (DM) halos, we assume that the population of 
these objects is well-described by a Press-Schechter distri- 
bution (e.g., Lacey and Cole 1993) 



dn(M, z) 
dM 



2 p(z) diner 
kTP dlnM 



v(M, z) exp 



\M,z) 



(1) 



where dn(M, z) is the number density of collapsed halos in 
the mass range [M, M + dM] and p(z) is the mean back- 
ground density at redshift z. The normalized fluctuation 
amplitude v(M,z) is defined as 5 c o(z)/a(M), where a(M) 
is the variance of the fluctuation spectrum filtered on mass 
scale M and 6 c o(z) is the linearly evolved overdensity of a 
perturbation that has collapsed and virialized at a redshift z 
(see Appendix A). We assume a scale-free power spectrum 
P(k) oc k n , so the variance can be written 



a(M) = <7i 5 M" 



(2) 
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where the subscript indicates that the normalization is to 
a mass of 10 15 h^M® (with h 50 = H o /50 km s _1 Mpc -1 ) 
and a = (n + 3) /6. Connection to the conventional normal- 
ization as within 8/i -1 Mpc spheres is straightforward 

ffw = [y( 16 ^o X ) 3 P<A>/10 15 h^] a a s = [1.19fi ] Q cr 8 (3) 

where p c = 3 H 2 /8ttG is the critical density. The power 
spectrum normalization deduced from cluster abundances 
follows the empirical fitting function 

<T8 = Cl^ 2 (4) 

where there is good, but not exact, agreement in the litera- 
ture on the values of Ci and C2 (White, Efstathiou & Frenk 
1993; Viana & Liddle 1996; Eke, Cole & Frenk 1996). For 
our calculations, we use c\ — 0.60, 

c 2 (fio) =0.36 + 0.31fio -0.28^0 (5) 
for open models with Ao = 0, and 

c 2 (fi ) = 0.59-0.16fi +0.06fio (6) 

for models with Qq + Ao = 1. These fitting functions were 
taken from Viana & Liddle (1996). There is a random un- 
certainy on <rg of approximately ±27% > * ne ma i n component 
of which results from the cosmic variance in the local clus- 
ter population. Since this is taken into account in a different 
manner later in the paper (section 3.3), we treat the above 
value as exact. 

We assume that the bolometric X-ray luminosity of clus- 
ters follows a power law in mass and redshift 

L x = L 15 M p (1 + z) s (7) 

over a range of 10 13 to 1O 16 M in mass and z < 2 in redshift. 
Here, and throughout the paper unless specified otherwise, 
the mass M is in units of 10 h^ Mq . Although assuming 
that this simple mass dependence holds over three orders 
of magnitude may seem unreasonable, the abundance of ob- 
jects drops quite sharply outside a central range of about 1.5 
decades for the flux limits that we are considering. High- 
mass objects are rare at any epoch, and low-mass objects 
quickly become invisible at larger redshifts; thus, any devi- 
ations from a power law outside this range will have little 
effect on our predictions. The intrinsic luminosity at fixed 
mass increases with redshift for s > 0. For reasons which will 
be made apparent, we focus attention on such "positive lu- 
minosity evolution" models, though it is straightforward to 
extrapolate our results to models with negative luminosity 
evolution. 

What values of the parameters p and s are expected? 
On dimensional grounds, we can write a scaling relation for 
the bolometric luminosity as 

L x oc J dVp 2 T 1/2 oc M 4/3 (l + z) 7/2 I(M, z) (8) 

where the virial theorem T oc M 2 ^ 3 (l + z) and the assump- 
tion of clusters as regions of fixed overdensity J dVp 2 ~ 
Mp(z) ~ M{1 + z) 3 with constant gas fraction is used to 
produce the scalings on the right hand side. Here I(M, z) is 
a form factor which retains the information on the mean 



internal density and temperature profiles of the clusters. 
Kaiser (1986) derived the above scaling under the assump- 
tion of self-similarity of the cluster population across both 
mass and epoch, so that I(M, z) = const. Although self- 
similarity may apply to the cluster popultion, it requires 
rather restrictive conditions — gravitational shock heating 
should be the dominant heating mechanism, cooling unim- 
portant, variations in gas fraction must be small, clusters 
must have similar internal structure, and so on. EH91 pre- 
sented empirical evidence against self-similar scaling in mass 
for Q = 1 models from the shape of the luminosity func- 
tion. (We return to this point below.) Kaiser (1991) and 
EH91 presented alternative models which invoked constant 
entropy either throughout the cluster gas or in the central 
core, respectively. The constant core entropy models of EH91 
yield p=ll/6 and s = ll/4. 

The approach taken here is to let the observations dic- 
tate appropriate values of p and s. The value of p then re- 
flects a composite mean description summarizing the mass 
sensitivity of cluster internal structure, gas fraction, cooling 
flows, efficiency of star formation and other microphysics. 
Variance in this relation is discussed below. The evolution 
parameter s encompasses time-dependent phenomena such 
as the changing overall density of the universe, the efficiency 
of radiative cooling, and heating of the ICM via gravitational 
collapse or supernova injection. Bower (1997) discusses the 
connection between cluster entropy and the evolution pa- 
rameter. The value s = 11/4 indicates that these processes 
are balanced; higher values indicate that heating mecha- 
nisms are dominant, while lower values indicate that cooling 
mechanisms are dominant. 

The above discussion pertains to the bolometric cluster 
luminosity. In practice, the X-ray luminosity within some 
range of photon energies, denoted by energy band E, is re- 
quired to connect with the observational data. The observed 
luminosity Lx{E) is a fraction of the bolometric 

L X (E)=L X f E [T(M,z),z] (9) 

where the factor /e [T(M, z), z] is found by numerically in- 
tegrating the Bremsstrahlung emissivity over the proper en- 
ergy range. As in EH91, we use an approximation to the 
Gaunt factor of g(E,kT) = 0.9(£7fcT) -0 ' 3 in this calcula- 
tion. The temperature of the cluster is related to the mass 
according to the equation fcT(keV) = 3.96M 2/3 (1 + z), a 
scaling law that is well supported by three-dimensional hy- 
drodynamic simulations (Evrard 1990a, b; Evrard, Metzler & 
Navarro 1996).^] Since the cluster temperature and emitted 
photon energy both scale as (1 + z) the redshift dependence 
for a fixed received energy band drops out and the X-ray 
luminosity fraction can be easily tabulated as a function of 
M or T alone. We use an energy range of 0.1 — 2.4 keV when 
comparing our model to the local luminosity function, and 
0.5 — 2.0 keV when predicting cluster number counts. 



* Although using this equation assumes that the cluster gas is 
fully virialized, the value of the band fraction is not strongly de- 
pendent on temperature and minor deviations from equilibrium 
will not significantly affect the results of this paper. 
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Figure 1. The BCS XLF of Ebeling et al.(1997). A change in 
the slope of the abundance function is evident around an X— ray 
luminosity of 3 X 10 44 erg s _1 . The solid line is the least-squares 
fit to a Press-Schechter form as described in the text. The dotted 
lines are the curves that correspond to L15 and p values at the 
vertices of the one-sigma error ellipse derived from the covariance 
matrix of the fit. The BCS XLF extends over a larger luminosity 
range than is shown here, but the additional data points do not 
help to constrain the parameters further. 

2.2 Constraints from the local luminosity 

The free parameters of our mass- luminosity model are L15, p 
and s. Since the extent of redshift evolution is still in ques- 
tion, we allow s to vary in the range 0.0 < s < 6.0. L15 
and p can then be determined by inverting equation ([?]) and 
expressing the Press-Schechter abundance as n(Lx , z)dLx ■ 
The resulting function is then fit to the BCS XLF of Ebel- 
ing et al.(1997), using an average redshift of 0.1 for the en- 
tire sample. Figure [jj shows the results of fitting a standard 
CDM-like model (Q—l, n = — 1, Ao = 0, and s = 3.5) to these 
data, with resulting least-squares fit values of p = 3.38±0.17 
and L15 = 2.50 ± 0.16 x 10 44 h~ 2 erg s _1 . 

Values of L15 and p for any given cosmological model 
can be similarly determined to roughly 10% and 5% accu- 
racy, respectively. Figure ^| displays the variation in Li 5 and 
p with Qo and n at a constant value of s = 3.5. The level of 
accuracy displayed is representative of all models. We ana- 
lyze each model individually and employ its best-fit values 
for these parameters in subsequent calculations. All models 
provide good fits to the BCS XLF, with reduced chi-squared 
parameters of (x 2 ) ~ 1.2 over seven degrees of freedom. 

The behavior of p in Figure 2 can be understood by 
examining the Press-Schechter abundance function, equa- 
tion (1), and the variation of as with Qo in open and flat 
models, equations (5) and (6). For a given set of cosmological 
parameters the form of the underlying mass function is corn- 



Figure 2. The best-fit values of L15 and p for open cosmological 
models, assuming an evolution parameter s = 3.5. The error bars 
correspond to the formal 1-cr uncertainty given by the diagonal 
elements of the covariance matrix. The overall shape of the curves 
and the level of uncertainty are similar in models which have 
different s or include a cosmological constant. The dotted lines are 
the limits on p that result if the observed luminosity-temperature 
relation is adopted as an additional constraint; this topic will be 
discussed in detail in §4.1. The numerical value of L15, defined as 
the luminosity at mass 10 15 /i^Mq, scales as hgQ . 



pletely fixed. For masses smaller than M* = (y/2ais /S c o) 1 , 
equation (1) behaves like a power law: n(M) ~ M~ 2+a . If 
M > M*, however, we have n(M) ~ exp[-(M/M*) 2a ]. As n 
increases, the exponential decay gets stronger and the power 
law gets more shallow, creating a stronger bend in the mass 
function at the transition region. In making the transforma- 
tion from n(M, 2) to n(Lx,z), we find that M ~ L\^ p , so 
a larger p will stretch out the luminosity function more and 
produce a shallower bend at Lx(M*). Indeed, this is what 
we see; larger values of n systematically require larger val- 
ues of p to reproduce the observed bend in the XLF. The 
primary effect of increasing Oo is to decrease as and thus 
a(M), strengthening the exponential cutoff in the high-mass 
regime and again causing a sharper bend at M*; the result 
is the expected increase of p with Q. The implications of this 
parameter and the meaning of the limits drawn in Figure 2 
are discussed further in section 4.1. 

The behaviour of L15 follows from this analysis, but 
is more complex. Since the XLF has a discernable bend at 
about 2 x 10 44 erg sec -1 , the value of M t can be expected to 
have a strong influence on the derived normalization. Clearly 
a larger value of n will lead to a smaller value of M*, since it 
only appears in the exponent l/2a. On the surface it appears 
that increasing £Iq should also cause to drop, since it is 
proportional to CT15, but this also has the effect of raising the 
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Figure 3. The EMSS X-ray cluster distribution as reconstructed 
by Nichol et al.(1997) The solid line represents the same best- 
fit "standard CDM" model that was shown in Figure 1, and the 
dotted lines represent 1-ct variations. 



value of p. Since Lx(M±) oc Li$M+ p , the required value for 
L15 can be expected to display a minimum where the two 
effects are balanced, ff n is zero, the dependence of M* on 
erg is greatly weakened and the change in p dominates the 
behaviour of L15. 

There are dependencies on s and Ao as well, but these 
are slight compared to those described above. Increasing s 
decreases L15 according to the factor (l + z) s folded into the 
BCS fit, and does not affect p at all. Introducing a cosmo- 
logical constant decreases L15 and p for low density mod- 
els, mainly because these universes require a larger S c o and 
therefore have a smaller K'h. 

The "1-ct variations" shown in Figure 1 are the curves 
obtained by using the vertices of the covariance matrix error 
ellipse to modify L15 and p. The uncertainties depicted in 
Figure 2 are the formal 1-ct error bars given by the diagonal 
elements of the covariance matrix. Since the functional form 
of the XLF is nonlinear, these limits should be taken as 
a reasonable estimates of the variance rather than precise 
measurements of a confidence level. 

At this point we have six parameters that describe our 
models: Qo, Ao, n, s, L15, and p. The first four are arbitrary, 
and the last two are determined from local observations. 
This is all the information we need to construct flux-limited 
statistics of the underlying Press-Schechter abundance for 
a given cosmology. We can make an additional test of the 
method by comparing our BCS-fitted distribution to the 
EMSS cluster sample (Henry et al.1992), recently reanal- 
ysed by Nichol et al.(1997), which has the advantage of be- 
ing binned in redshift. The results of this comparison are 
presented in Figure ^, along with the same variations on the 



theoretical XLF that were displayed in Figure 1. The theory 
with s = 3.5 and data are in good agreement, but the depth 
of sample is too shallow to provide significant leverage on 
s; models with s — 0.0 or s = 6.0 provide equally good fits. 
Putting good constraints on cluster evolution requires either 
a very large sample or a very deep sample, as we will see in 
the next section. 



3 PREDICTIONS AND OBSERVATIONS 
3.1 The logN-logS 

The model as stated above is sufficient to describe the un- 
derlying abundances of clusters for a variety of universes. 
We wish to compare these abundances to recent observa- 
tions of the logN-logS statistic. Calculating number counts 
as a function of observed flux is done by integrating the 
number density of clusters over mass and redshift, with a 
lower limit on the mass integral coming from the minimum 
mass required to meet the flux level at each redshift. 

The band- limited flux Sx (E) is related to the observed, 
band limited X-ray luminosity by the usual relation 



S X (E) 



Lx(E) 
47rr 2 (l + z) 



(10) 



where r is the physical distance from the observer to the 
cluster determined by the cosmological parameters within 
the Robertson-Walker metric. For a universe with Ao = 0, 
this gives 



Sx(E) = 



L 15 H$qS(l + z) a M p f E (M) 
47rc 2 [ 9 „2 + (g„ - l)(V2q 2 + l-l)] 2 



(11) 



This equation allows the minimum mass satisfying a given 
flux limit to be calculated at each redshift. 

Depending on the form of the observations that we are 
trying to imitate, we have the option of adding in at this 
point a correction for the finite cluster size and the point 
response function of the telescope, or some other factor that 
represents the efficiency of the flux recovery. This factor can 
either be a pure number or, if you are willing to invoke a 
model for the surface brightness, a function of the cluster's 
angular diameter (e.g., EH91). As the observers have al- 
ready made this correction in constructing the logN-logS, 
we calculate the total flux as given above. 

The number of clusters per unit mass on the sky is 
calculated from the basic relation 



dN = dn{M, z)dV = dn(M, z) 



r dr 



VI - kr 2 



dQ 



(12) 



which can, with the appropriate function r(z), be trans- 
formed into a function of redshift 



(13) 



dN _ dn(M,z) [q z + (go -l)(y/2q z + 1-1)] 
dzdQ H$(l + z) 3 q% ^1-250 + 290(1 + 2) 

which again holds for Ao = 0. Integrating this function over 
redshift with lower mass limit given by equation ([ll]) gives 
the number of visible clusters in the sky per steradian as a 
function of limiting flux Sx{E). Similarly, we can construct 
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a differential logN-logS or keep track of the total number in 
each redshift bin to predict redshift distributions. 

The situation for fio + Ao = 1 universes is slightly more 
complicated in that the function r(z) can not be expressed 
in simple analytic form. It can, however, be tabulated and 
used to calculate the volume element and observed flux to 
arbitrary accuracy. Appendix B contains the details of this 
derivation. 

3.2 ROSAT Data and Comparisons 

We use logN-logS data from two independent, serindipitous 
samples of X-ray clusters derived from deep, pointed ROSAT 
observations. The work of Rosati & Delia Ceca (1997) in- 
cludes 125 clusters over a total sky area of about 35 square 
degrees, and goes down to a flux limit of 2 x 1CP 14 ergs sec 
cm' 2 . The work of Jones et al.(1997) includes 34 clusters 
and goes from 4 x 10" 14 to 2 x 10" 13 in Sx- Both groups 
employ different methods and assumptions which we now 
discuss. 

Rosati & Delia Ceca (1997) use a wavelet decomposi- 
tion algorithm to identify extended sources in ROSAT PSPC 
fields, and employ the wavelet coefficients to reconstruct the 
total flux of each source. This technique is subject to large 
random errors on clusters of low signal-to-noise ratio, but is 
no worse off in these terms than fitting the cluster profile 
to a standard (e.g., King) model. There is also a bias inher- 
ent in reconstructing the area under a non-gaussian profile 
with a gaussian wavelet, but this is small compared to the 
sources of random error in the problem. In his doctoral the- 
sis, Rosati (1995) presents details of the reconstruction pro- 
cess. He also implements a sophisticated correction for the 
sky coverage, which takes into account distortion in PSPC 
images at large off-axis angles and varies with the angular 
diameter and flux of the source. Uncertainties in the inte- 
gral counts include Poisson noise, uncertainties in the sky 
coverage fraction, and the random errors inherent in recon- 
structing the flux of an image with low signal-to-noise ratio; 
the error bars used in the differential counts, however, are 
just Poisson. 

The work of Jones et al.is derived from the WARPS 
cluster survey of Scharf et al.(1997), which uses the VTP 
analysis method (Ebeling & Wiedenmann 1993; Ebeling 
1993) to detect sources. To account for any flux that might 
be hidden under the background level, they assume a stan- 
dard /3-model form with (3 = 2/3, derive a normalization 
and core radius from the source profile, and integrate in 
the region of low signal-to-noise ratio to find an appropriate 
correction factor. The model is not used to determine the 
overall flux, just to make a second-order correction to the 
VTP count rates. Details of their flux correction procedure 
can be found in Ebeling et al.(1996). The error bars of this 
dataset are strictly Poisson and are claimed to dominate all 
other sources of error. Numbers used in this analysis were 
kindly made available by L. Jones prior to publication. 

Both groups attempt to construct the logN-logS in a 
largely model independent fashion. The advantage of this 
approach is that it avoids building in biases based on ques- 
tionable assumptions. The disadvantage is that it makes it 




-2 -1 



log S (0.5-2.0 keV) erg sec-' cm- 2 ] 

Figure 4. The integrated cluster surface density above an X- 
ray flux limit S. The solid line respresents the best-fit values of 
Lis = 3.556 and p = 2.174, while the dashed lines represent 1- 
a variations around these values as described in the text and 
displayed in Figures 1 and 3. 

difficult to assess sample completeness in a systematic fash- 
ion; it is unclear how large a population may be missing 
from the detected counts because of low surface brightness 
and/or large angular extent. Of course, comparing results 
of the independent groups is a simple gauge of systematic 
effects. The data, shown in Figure 4, indicate consistency 
within modest statistical errors. Given that it is easier to 
err in the direction of missing extended, low surface bright- 
ness sources, it is reasonable to interpret the current data 
as providing firm lower limits to the true cluster counts, 
and accurate estimates if no such population exists. Thus, 
a model which underpredicts the number of clusters can be 
excluded with somewhat greater confidence than one which 
overpredicts the observed counts. 

Figure ^ gives an example of the integrated logN-logS 
for a universe with Qo — 1, Ao = 0, s = 3.5, and n = —2. Su- 
perimposed are the observational data from Rosati & Delia 
Ceca (1997) and the WARPS sample. The data are shown 
in their familiar cumulative form here, so the points of each 
sample are not statistically independent from one another. 
A crude judgement of goodness of fit can be made by com- 
paring the model to the faintest points of each sample. This 
particular model underpredicts the number of clusters, to a 
moderate but still significant degree. This is true even for 
the logN-logS curves that arise from modifying L15 and p 
through the vertices of their one sigma error ellipse, shown 
as the dashed lines in the figure. 

Figures 5 through 7 demonstrate the sensitivity of the 
counts to variations in Qo and s. Figure [| displays the 
changes that occur as Qo decreases from 1.0 to 0.2 while 
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Figure 5. The integrated cluster surface densities for models with 
Ao = 0, n = — 2, and s = 3.5, while Qo takes on the values 0.2, 0.4, 
0.6, 0.8, and 1.0 



Figure 6. The integrated cluster surface densities for models with 
Ao = 0, Ho = 1.0, and n = — 2, while s takes on the values 0.5, 1.5, 
2.5, 3.5, 4.5, and 5.5. 



s is held constant at a value of 3.5. Clearly, low f2 universes 
produce higher counts, a reflection of the slower evolution 
experienced by the cluster population in such models (Rich- 
stone, Loeb & Turner 1992) as well as the larger volume 
element per redshift interval. Figure Q displays the logN- 
logS for the same five universes but with a mild evolution 
parameter s= 1.0. Reducing s has the effect of reducing the 
differences between universes with different density parame- 
ters as well as lowering the overall number of cluster, since a 
lack of positive luminosity evolution quickly dims the high- 
redshift population expected in open models. This extinction 
has a proportionally larger effect on the low fio models, for 
which a given flux limit represents a deeper probe. 

The effect of changing s while holding flo fixed is shown 
in Figure [| The variation in counts seen here is a combina- 
tion of two factors: first, a model with less positive evolution 
in the cluster luminosity will contribute fewer objects to the 
logN-logS at a given flux limit; and second, the overall nor- 
malization of the mass-luminosity relation will be larger if s 
is smaller. (The second effect arises from the approximation 
that the entire BCS sample lies at an average redshift of 0.1, 
which introduces a dependence of L15 and p on s-see equa- 
tion (^).) The two effects change in strength and work in 
opposite directions, so the variation shown in Figure 6 can- 
not be taken as universal. In particular, low-density models 
display a much stronger dependence on s than models in 
which £7o = 1; this can be understood as an increase in the 
importance of the high-redshift population relative to the 
local sample. 




-2 -1 

log S (0.5-2.0 keV) [1Q-' Z erg sec' cm Z J 



Figure 7. The integrated cluster surface densities for models with 
Ao = 0, n = — 2, and s = 1.0, while Qo takes on the values 1.0, 0.8, 
0.6, 0.4, and 0.2. Note that a lower value of s greatly supresses the 
differences between low-density and high-density models, because 
high-redshift clusters become harder to observe. 
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Figure 8. The square root of the reduced chi-squarcd parameter 
for all models is displayed, as well as the "uncertainty" in this 
parameter, the derivation of which is described in the text. In the 
upper row, dottcd-line contours have been drawn corresponding 
to (x 2 ) 1 ^ 2 = 1-07, 1.28, and 1.51, which represent formal proba- 
bilities of about 30%, 5%, and 0.3% that the data come from the 
same underlying distribution as the model and differ only through 
gaussian random variations. Solid-line contours have been drawn 
at (x 2 ) 1 ^ 2 = 2.0, 3.0, and 4.0, which correspond to very small 
probabilities but may be useful in combination with the uncer- 
tainties. The dashed lines are drawn every time the likelihood 
parameter doubles thereafter, i.e. at levels of 8.0, 16.0, etc. In 
the lower row, the dotted-line contours are drawn at uncertainty 
levels of 0.1, 0.2, and 0.3; the solid lines at levels of 0.4, 0.6, 0.8, 
and 1.0, and the dashed lines at each doubling over 1.0. Note 
that although the uncertainties get very high in the region of low 
Qo and high s, the likelihood parameter is larger still by about a 
factor of 4. Figure (8a) displays models with n = 0, Figure (8b) 
with n = — 1 and Figure (8c) with n = — 2. 

3.3 Examination of the Qq — s plane 

Changing the cosmological parameters over the proposed 
range can produce massive overabundances or moderate un- 
derabundances. To assess the likelihood of each model we 
compare the predicted differential logN-logS (the number 
of clusters in a specified flux bin) to the data points and cal- 
culate the reduced chi-squared factor for that model. The 
full set of data used to constrain our parameter space in- 
cludes the faintest 6 data points from the WARPS sample 
and ten points from Rosati's, for sixteen degrees of free- 
dom altogether. The use of differential data means that sta- 
tistical errors are formally independent. We employ a chi- 
squared statistic because it provides a measure of the ab- 
solute goodness of fit for the models, assuming the error 
bars are Gaussian. Our results would be similar were we to 
employ a maximum likelihood approach. The main disad- 
vantage of this approach is the difficulty in including the 




X = Q + X= 1 




uncertainty in L15 and p for a particular model. We adopt 
an approximate procedure to solve this problem, which will 
be described shortly. 

Figures 8(a-c) present the "likelihood" plots of all mod- 
els having n = 0, — 1 and —2, respectively. The value dis- 
played in the upper row for each model is the square root 
of the reduced chi-squared parameter (x 2 ) 1 ^ 2 > where the av- 
erage is over all data points in the logN-logS. The dotted 
contours in the upper row mark models which deviate from 
the data at one, two, and three-sigma confidence levels. The 
main limitation in the ability of theory to predict abun- 
dances accurately at this level of analysis is the determina- 
tion of Lis and p. Although an accuracy of 5% to 10% is 
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quite good (and quite possibly as good as we're going to 
get, since this method is limited by the number of nearby 
clusters) , varying these numbers within their confidence lim- 
its can still produce considerable changes in the predicted 
logN-logS (see Figure 4). 

To account for this, we calculated the likelihood pa- 
rameter, (x 2 ) 1 ^ 2 , f° r the four vertices of the standard error 
ellipse in L 15 -p space as well as for the best-fit model. The 
uncertainty shown is the difference in (x 2 ) 1 ^ 2 between the 
best-fit case and the perturbed case which gave the small- 
est value of the likelihood parameter. The results are pre- 
sented in the lower rows of Figures 8(a-c) and can be inter- 
preted as a "one-sigma" uncertainty in the contours of the 
top row. Keep in mind that this variance is not symmet- 
ric; we present only the difference in the direction of greater 
likelihood. Since our aim here is to provide a reasonable in- 
dication of the uncertainty in the modeling, this heuristic 
error estimate seems sufficient. The data were presented in 
this manner to make it possible for the reader to estimate 
constraints at other confidence levels. 

ft is clear from Figures 8(a-c) that the results are not 
strongly affected by the value of Ao . Introducing a cosmolog- 
ical constant increases the volume element per redshift bin 
throughout most of the integration range, which reduces the 
observed flux of each cluster, but it also increases the total 
surface density of clusters. As can be seen from a close ex- 
amination of the figures, in most regions of the parameter 
space introducing a cosmological constant slightly raises the 
value of fio required for a good fit. Since raising fio lowers 
the number of clusters, we conclude that increased integra- 
tion volume is the dominant effect for most models. In other 
words, setting Ao 7^ tends to increase the surface density of 
clusters, and a higher value of fio is needed to compensate. 

Changing the spectral index n has a quite strong ef- 
fect on the derived mass- luminosity relationship (Figure 2), 
which as we will see can have important consequences. Be- 
cause all of the models are constrained to agree with local 
abundances, however, n has limited leverage on the distant 
cluster population. The main effect on the logN-logS is that 
models with larger n have sharper exponential cutoffs in 
their mass distribution, resulting in a lower cluster surface 
density. The best-fitting models for n = thus have lower fio 
and greater luminosity evolution to compensate. 

The clearest relationship visible in Figures ^|(a-c) is, of 
course, that between the density and evolution parameters, 
fio and s. It seems from these results that if we wish to 
believe in the degree of evolution indicated by self-similar 
or constant entropy arguments — s = 3.5 and 11/4, respec- 
tively — then we obtain a useful lower bound of fio > 0.3. 
The lower end of this range is consistent with limits from 
the mean intracluster gas fraction of clusters, which give 
fi /i 2/3 = 0.28 ± 0.07 (Evrard 1997) in the most straightfor- 
ward interpretation. 

Conversely, if the degree of luminosity evolution is in- 
stead very small — s < 1 — then the analysis strongly rules 
out a critical density universe. Recall that all the off-ridge 
models which have fio = 1 underpredict the number of ob- 
served clusters (see Figures 5-7), and these should be treated 
more harshly given the possibility for survey incompleteness. 



It seems that if the universe does indeed have a critical den- 
sity, we need at least a moderate degree of evolution in the 
cluster luminosity to fit the data well. The assumption of 
constant cluster entropy (s = 11/4) put forward in EH91 
would be sufficient to allow such models. The preferred re- 
gion, however, indicates an even larger degree of evolution. 
This could come about if radiative cooling occurs on a much 
longer time scale than gravitational collapse (Bower 1997), 
or if a significant amount of energy is injected into the ICM 
by galaxies (e.g., Metzler & Evrard 1994). 

The results of this analysis are summed up in Table |l[ 
The limits put down are at the 95% (about 2 a) confidence 
level (i. e. there is only a 5% chance of the model randomly 
producing the observed data), after reducing the (x 2 ) 1 ^ 2 pa- 
rameter by the uncertainty indicated in Figures 8(a-c). The 
models are examined in intervals of 0.05 in fio an d 0.5 in 
s, so a > or < sign just means that the limiting model 
was at or very near the required confidence level of 95%. 
In addition, all of the disallowed models indicated here un- 
derpredicted the number of clusters, except for those which 
assumed s = 3.5. The constraints thus allow for a one-sigma 
variation in L15 and p tailored to bring the model closer to 
the data, and can be treated as conservative conclusions. 



4 ADDITIONAL CONSTRAINTS FROM THE 
CLUSTER POPULATION 

4.1 The Luminosity- Temperature Relation 

Because of the degeneracy between intrinsic luminosity evo- 
lution and cosmological evolution, the logN-logS alone limits 
models to the ridge seen in Figure ^[ It is worthwhile to see 
what additional constraints can be placed on the parame- 
ter space by including additional observational information. 
One interesting question is what each model predicts for the 
X-ray luminosity — temperature correlation. So far, we have 
made use of the virial relationship T oc A/ 2//3 only to tab- 
ulate the band fraction. Because the fraction fs{T) for the 
ROSAT bands is a weakly dependent function of tempera- 
ture, our results up to this point would change very little 
if we were to adopt a different T(M) behavior. If we are 
willing to promote virial equilibrium to a strong assump- 
tion, however, we can use it as an independent test of the 
parameter p. There is good theoretical support from numer- 
ical simulations for virial equilibrium within the non-linear 
portions of clusters (Evrard et al.1996) and modest empir- 
ical support for this assumption from analysis of the mean 
intracluster gas fraction (Evrard 1997). Comparison with 
the values derived from the BCS XLF fit then provides an 
added, non-trivial constraint on the models. 

The virial scaling assumption T oc M 2//3 (l + 2) trans- 
lates directly into a bolometric Lx~T relation of the form 

L x oc T 3p/2^ 1 + z y-3 P /2 ^ 

The observed relationship is, roug hly, L x oc T 2 8±0 - 3 . (Edge 
& Stewart 1991; Arnaud 1994). The quoted error in the slope 
is more generous than those reported in individual works, in 
order to allow for possible systematic uncertainties betwen 
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s 
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s = 3.5 


n > 0.20 


fi > 0.25 


H > 0.35 


n > 0.20 


C > 0.25 


n 


> 0.35 


s = 1.0 


n < 0.45 


Q < 0.60 


f^o < 0.75 


C < 0.45 


tt < 0.65 


n 


< 0.80 


s = 0.0 


n < 0.30 


n < o.4o 


H < 0.60 


n < 0.35 


n < 0.45 


n 


< 0.65 



Table 1. A summary of the constraints that can be placed on f!o an d s using the information in Figures 8. All of these constraints are 
at the 95% confidence level, after accounting for the displayed uncertainty in the chi-squared parameter. 



different data sets. Consistency with the observed slope re- 
quires p = 1.9 ±0.2. This range, drawn in Figure 2 as dotted 
horizontal lines, incorporates the value 11/6 appropriate for 
the constant central entropy model of EH91 (see also Bower 
1997). A quick glance at the figure reveals that including this 
result puts a strong constraint on f2o and n. If n — — 1, this 
range of values is found only in very low-density universes 
(Q,o < 0.1). For models with n = —2, on the other hand, we 
find acceptable values of p in the range 0.25 < flo < 1.0. If 
a cosmological constant is included, the allowed range in Qo 
shifts to slightly higher values. 

The real perturbation spectrum is not likely a pure 
power law, but this analysis indicates models with effective 
spectral index n <C —1.5 are favored in universes with rea- 
sonable values of the density parameter. Such "red" spec- 
tral values are also favored by the shape of the tempera- 
ture abundance function (Henry & Arnaud 1991; Oukbir, 
Bartlett & Blanchard 1997). For reference, Huss, Jain & 
Steinmetz (1997) present a list of effective spectral indices 
(n e // = dln P(fc)/dln k) for a number of popular cosmolog- 
ical models. A cold, dark matter model with fio = 0.3, A = 
and h = 0.7 has n e // = — 1.74, for example. 

Room to maneuver around this conclusion can be 
gained by invoking a degree of scatter in the relationship be- 
tween bolometric X-ray luminosity and binding mass. Scat- 
ter will tend to flatten the slope and increase abundance at 
the bright end of the predicted luminosity function, implying 
a bias toward overestimating p and Lis compared to the case 
with no scatter. The amount of scatter is not known a priori, 
but a likely lower bound can be obtained from gas dynamic 
simulations. The two sets of 18 cluster simulations of Metzler 
(1995) revealed a scatter of magnitude (<5(log Lx) 2 ) = 0.047 
for simulations involving dark matter and baryon fluids, and 
(<5(log Lx) 2 ) = 0.026 for simulations which also included a 
model for the ejection of gas from galaxies into the ICM. 
We have checked that including scatter at the larger of the 
above values results in a fractional change in the best fit p of 
at most 10%, comparable to the uncertainty from the BCS 
XLF fit. We conclude that such modest levels of scatter are 
not strongly affecting our results. Larger amounts of scatter 
would lead to downward revisions in our best fit values of 
p in Figure 2, and so would somewhat relax constraints on 
the effective spectral index on cluster scales. 

Another important prediction from equation ( |l4| ) is that 
we expect redshift evolution in the bolometric luminosity at 
a fixed temperature proportional to (1 + z) s ~ 3p//2 . Since fit- 
ting the local relation requires 3p/2 — 2.8 ± 0.3, the expec- 



tation is of an offset in the intercept of the Lx — T relation 
of amplitude 

A (log L x (z))\t = [s- (2.8 ±0.3)] A(log(l + *)) (15) 

independent of spectral index or cosmology. Both the orig- 
inal self-similar value s = 7/2 and the constant central en- 
tropy value s = 11/4 predict very little change in the inter- 
cept of the Lx~T relation back to z~0.5. The full range we 
investigate, s £ [0, 6], translates into shifts of factors ranging 
from about one-third to three at z = 0.5. 

Comparing a recent compilation of 15 distant clus- 
ters ({z) ~ 0.33) by Henry et al.(1997) with a similar 
sample of the nearby population ({z) ~ 0.07) by David 
et al.(1993) provides us with an early estimate of this 
shift. Comparing the intercepts of the two samples, we find 
A(log Lx)/A(log(l ± z)) = 1.1 ± 1.1 at the one-sigma level, 
which corresponds to a plausible range of 2.5 < s < 5.3. 
The observations rule out very low values of the evolution- 
ary parameter s, and therefore eliminate the low Qo end of 
the allowed ridge in Figure H. 



4.2 Redshift Distributions 

Additional constraints can be found by examining the red- 
shift distribution of the clusters, which can be expected to 
vary strongly with our choice of cosmology. Indeed, Ouk- 
bir & Blanchard (1997), employing methods similar to ours, 
found the redshift distribution in the EMSS sample (Gioia 
& Luppino 1994) to be well fit by either a low density model 
with mild evolution (f2o=0.2, s = 0.5 in our terminology) 
or an Einstein-deSitter model with more rapid evoluition 
(Oo = 1, s = 3.8). These models lie along the ridge allowed 
by the ROSAT counts in our analysis. 

As discussed in section 3.1, the redshift distribution 
is straightforward to calculate for a given model and flux 
limit. Although obtaining accurate redshifts for clusters is 
a time-consuming process, the shape of this statistic is not 
strongly dependent on small differences in the efficiency of 
flux reconstruction or the choice of model for cluster sur- 
face brightness. A confident flux limit and a well-understood 
sky coverage function, are, however, necessary to predict the 
normalization of these curves. 

Figure ^| displays the differential and cumulative dis- 
tributions at a flux limit of 5.5 x 10~ 14 erg sec -1 cm~ 2 in 
the band 0.5-2.0 keV, the parameters of the WARPS sur- 
vey. Each category of models is represented by four specific 
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Figure 9. The cumulative and differential rcdshift distributions 
for the A = 0, n = —2 models described in the text. The effect of 
changing Qo and s can be clearly seen here; increasing Qo sharply 
reduces the number of high-redshift clusters, while increasing s 
brings the lower-mass objects above the flux limit. 

choices within the region of plausible likelihood, all of which 
assume n = — 2: 

(i) Oo = 0.3, s = 1.0 

(ii) Q = 0.3, s = 2.75 
(hi) Qo = 1.0, s = 2.75 
(iv) Q = 1.0, s = 6.0 

Models (i) and (iv) fit all of the logN-logS data quite 
well, lying near the center of the ridge described by Figure 
8c, but model (i) predicts a decrease in the zero-point of 
the L x -T relation larger than observed. Models (ii) and (iii) 
are given as reference points, and to display the effects of 
moving perpendicular to the ridge in the Slos plane. These 
models are ruled out at about the 95% confidence level after 
subtracting the uncertainty in their likelihood parameters, 
but make the effect of changing the parameters clearer. 

As you can see, the redshift distribution does not dis- 
criminate well between models which lie on the ridge in Fig- 
ure 8c, but displays significant differences when the param- 
eters are changed in other directions. Samples of size ap- 
proaching 100 clusters are soon arriving (Rosati 1997, pri- 
vate communication) and these will provide stringent tests of 
the models by constraining the area under the high redshift 
tail of the distribution. In particular, if Qo — 0.3 as indi- 
cated by the mean intracluster gas fraction and s is close to 
the 11/4 fixed core entropy value (model ii), then we predict 
10% of clusters above the WARPS survey flux limit should 
lie at z > 0.75. 



5 CONCLUSIONS AND FUTURE 
DIRECTIONS 

We assume that the underlying number density of clusters as 
a function of mass is well described by the Press-Schechter 
function and model the scaling of cluster bolometric lumi- 
nosity and temperature as power-laws in binding mass and 
epoch. For a given cosmological model, we determine the free 
parameters of the current mass-luminosity relation by fitting 
the predicted X-ray luminosity function to the BCS XLF 
(Ebeling et al. 1996), and integrate the abundance function 
out to high redshifts to compute counts. 

A direct comparison by means of chi-squared analysis 
made to differential logN-logS data from two ROSAT sur- 
veys produces a ridge of acceptable models following roughly 
s ~ 6fio. If s is in fact about equal to 3 as expected from 
simple theoretical arguments (Kaiser 1986, EH91) then the 
universe seems to require at least a moderate density param- 
eter Oo ^ 0.3 in order to bring the cluster number counts 
down to observed levels. Conversely, a universe with crit- 
ical density is inconsistent with a low evolution parameter 
s < 2.5; such models are strongly ruled out because they un- 
derpredict the number of clusters and are harder to save by 
invoking systematic errors in the surveys. Although chang- 
ing the slope of the power spectrum and incorporating a 
cosmological constant also affected these results, these pa- 
rameters contribute only slightly to the overall behaviour. 
An independent constraint from the local, observed Lx~T 
relation limits singles out the n = — 2 plane, further reducing 
the space of allowed models. 

The logN-logS alone is a useful, but limited, discrimi- 
nator of cosmological models. Because the local XLF is used 
to determine the mass-luminosity relationship, the predicted 
counts will always be limited in accuracy by the level of shot 
noise in the local cluster population. Assuming that the lo- 
cal XLF will not increase much in accuracy, this method has 
the potential to pin down a region in the Qo-s plane with a 
width of about ±0.1 in fio and ±1.0 in s. This estimate was 
made under the assumption that the uncertainties arising 
from the XLF are the only significant ones; in practice this 
level of discriminating power could be achieved when the 
overall uncertainties on the logN-logS data have decreased 
to about one-third their current size. 

One could try to improve the accuracy of this method 
by pushing observations down to a lower flux limit, but the 
returns are limited by the local XLF uncertainty. In order 
for a new data point at a fainter flux limit to be useful, it 
needs to have error bars significantly smaller than the un- 
certainty in the logN-logS from the local XLF. For example, 
the uncertainty on the integrated logN-logS at a flux limit of 
10~ 15 erg s" 1 is about 30%. We estimate that a data point 
with 10% error bars would make a useful contribution to 
the discriminating power of this method. Assuming simple 
Poisson errors, this level of precision could be achieved with 
a sample of about 100 clusters, over 1-2 square degrees of 
sky. Toy models show that such a point would strengthen 
the lower limit on s for critical density universes, and greatly 
increase the lower limits on Qo for models with moderate s. 

Independent observations can be used to narrow the al- 
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lowed range in parameter space. Assuming clusters are in 
virial equilibrium, the shift in the zero-point of the L x -T 
relation limits the evolutionary parameter s. Current obser- 
vations at (2} ~ 0.3 (Henry et al.1997; Mushotzky & Scharf 
1997) indicate s lies in the range 2.5 to 5. Similar data 
at higher redshifts, with subsequently longer lever arm in 
equation (|l^), would be very useful in reducing the allowed 
range. Redshift distributions of X-ray flux limited samples 
provide additional independent constraints (Oukbir & Blan- 
chard 1997; Bower 1997). The extent of the high redshift tail 
of the distribution is a sensitive indicator of which side of the 
Q.Q—S ridge the universe lies. If Qo — 0.3 as indicated by the 
mean intracluster gas fraction (Evrard 1997) and s > 11/4, 
then we predict at least 10% of clusters above the WARPS 
survey flux limit should lie at redshifts z > 0.75. 
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APPENDIX A: CALCULATING 5 C {Z) 

The critical overdensity S c (z) is defined as the linearly ex- 
trapolated overdensity of a perturbation that has just col- 
lapsed at redshift z. For universes with Oo = 1, this num- 
ber is the canonical 8 C — 1.686; for open universes or flat 
universes with Ao 7^ this overdensity will be somewhat 
smaller. Since the CDM power spectrum used in our analy- 
sis is normalized to the COBE microwave background mea- 
surements and therefore the power spectrum at the present 
day, the critical overdensity of a perturbation collapsing at 
redshift z is further extrapolated to a redshift of zero us- 
ing the linear growth factor appropriate to the cosmology. 
This quantity is labeled S c o(z). All the following relations 
are derived using the spherical collapse model. 

For a flat universe with Qo = 1, we use the following 
relations 
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S c (z) 
5 c o{z) 



-(12^) 2/3 
20 V ' 

{l + z)S c (z) 



(Al) 
(A2) 



For an open universe with flo < 1 and Ao = we use 
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the relationship derived by Cole and Lacey (1993) 

2/3 



D(0) 
D(z) 



Sc(z) 



5 c (z) = 
Sco(z) = 
D(0) = 

"0 Xq' ~ 

where xo = fig 1 — 1,1] = cosh -1 (2/Q,{z) — 1), and D repre- 
sents the linear growth factor. 

Finally, for a flat universe with flo + Ao = 1 we use 
an approximate parameterization to the Eke et al. (1996) 
results for S c (z) 



1 _3_ 3 VI + xo 

xo 



(A3) 
(A4) 

, ln(Vl + x - V^) (A5) 



6 c (z) = 1.68660(1 +0.01256 log Q(z)] 



(A6) 



The form of the parameterization was inspired by Kitayama 
and Suto (1997). This reference also contains the functional 
form of the exact solution, which is a hypergeometric func- 
tion of type (2,1). To get & c o{z) for this type of universe we 
use the solution for the linear growth factor found in Peebles 
(1980), §13 



r a .3/ 2(a .s + 2) -s/ 2da . i 
x Jo 



(A7) 



where x = a/a e , and a e = [(1 — Ao)/(2Ao)] 1 ^ 3 , the inflection 
point in the scale factor. In our analysis this function is 
integrated numerically to find the growth factor at redshifts 
z and 0, and then S c o(z) again equals 5 c (z)D(0)/D(z). 



solution in terms of elliptic integrals of the first kind, but 
expressing it in this manner requires the integral to be per- 
formed from a = 0. In addition, for Ao greater than about 
0.7, the elliptic integral must be evaluated outside of its de- 
fined domain. It is computationally more reliable to evaluate 
the integral numerically. For those who are interested, the 
"exact" solution is 



r(z) 
ro 



1 + (1 - \/3) 



Uo 



- F 



1 + (1 + s/S)uo 
1 + (1-V3)«z 



1 + (1 + V3)t 



where u = (A /fio) 1/3 , u z = u /(l + z), k = \J 2 + V3/2, 
and r = c/(H nl/ 3 \l /fi V3). 



APPENDIX B: THE VOLUME ELEMENT FOR 

Ao ^0 

The number of clusters in a given volume element dV cen- 
tered at redshift z and mass range dM is given by the Press- 
Schechter abundance times the differential volume 

r 2 rlrdQ 

dN(M) = n(M, z)dVdM = n(M, z) —= dM (Bl) 

VI — kr 2 

To integrate this function over a range in redshift rather than 
in physical distance, we need to find r(z) for a particular 
cosmology and make the transformation to the form 

dN(M) = n(M, z)f(z)dz dQ dM (B2) 

For a universe with Qo + Ao = 1, k = and f(z) = r 2 dr/dz. 
The function r(z) can be tabulated by integrating along the 
null geodesic: 

(B3) 



r , n da 

r = / dr = — 

Jo Ja a(1 



We can then express the Friedmann equation in terms of 
today's values of the cosmological parameters 



d = a# y^ + Ao (B4) 

and insert this into the geodesic integral to get r(z) in a 
straightforward manner. The integral obtained has an exact 
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